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Abstract: The search for Galactic pevatrons is now a well-identified key science project of all instruments 
operating in the very-high-energy domain. Indeed, in this energy range, the detection of gamma rays 
clearly indicates that efficient particle acceleration is taking place, and observations can thus help identify 
which astrophysical sources can energize particles up to the ~PeV range, thus being pevatrons. In the 
search for the origin of Galactic cosmic rays (CRs), the PeV range is an important milestone, since the 
sources of Galactic CRs are expected to accelerate PeV particles. This is how the central scientific goal 
that is ‘solving the mystery of the origin of CRs’ has often been distorted into ‘finding (a) pevatron(s)’. 
Since supernova remnants (SNRs) are often cited as the most likely candidates for the origin of CRs, 
‘finding (a) pevatron(s)’ has often become ‘confirming that SNRs are pevatrons’. Pleasingly, the first 
detection(s) of pevatron(s) were not associated to SNRs. Moreover, all clearly detected SNRs have yet 
revealed to not be pevatrons, and the detection from VHE gamma rays from regions unassociated with 
SNRs, are reminding us that other astrophysical sites might well be pevatrons. This short review aims at 
highlighting a few important results on the search for Galactic pevatrons. 
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1. Introduction 
1.1. The Hunt for the Sources of Galactic Cosmic Rays 


The term Pevatron, simply referring to an object capable of accelerating particles up 
to the PeV (=10' eV) range, is now widely used in the scientific programs of all major 
instruments investigating particle acceleration in astrophysical sources. This passionate search 
for pevatrons is not a stubborn search for the most powerful astrophysical particle accelerators, 
but is undoubtedly motivated by the hunt for the origin of Galactic cosmic rays (CRs). 

CRs are relativistic charged particles, typically ~92% protons, ~6% Helium, ~1% elec- 
trons and ~1% heavier nuclei, that fill the whole disk of the Milky Way. Although decades 
of experiments and theoretical studies helped to accumulate valuable knowledge, the funda- 
mental question of their origin is still missing [1]. 

It is widely accepted that accelerators located inside our Galaxy must produce the bulk 
of CRs detected at the Earth. One strong argument for the Galactic origin comes from the 
detection of gamma rays: observations of the Galaxy revealed that the Galactic disk is gamma- 
ray bright, which is a clear indication of the interactions of CRs with interstallar medium (ISM) 
material and Galactic magnetic field in the disk, and that gamma rays roughly scale with the 
amount of material in the disk. Moreover, observations of the Large Magellanic Cloud (LMC), 
located in the vicinity of our Galaxy at ~50 kpc, have shown that the gamma-ray signal is 
weaker than expected given the mass of the target material in the LMC [2,3]. This reduced 
gamma-ray signal from the LMC can be seen as evidence that CRs are not homogeneously 
distributed in the Universe, and originate within our Galaxy. Additionally, a simple reasoning 
on the Larmor radius rz of CRs, compared to the typical Galactic halo H size imposes that 
CRs of energy E remain confined within the Galaxy if rz (E) < H. For typical magnetic fields 
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of the interstellar medium (ISM) of the order of ~uG and a halo size of a few ~kpc (although 
both of these quantities remain poorly constrained), this simple estimate indicates that CR 
protons of energy E < 10!°-!”eV are expected to remain confined in the Galaxy. Therefore, 
we have indications that (1) CRs originate inside our Galaxy; (2) CRs are confined inside our 
Galaxy up to at least < 10'°-10!eV. These CRs of Galactic origin are simply called Galactic 
CRs. 

Moreover, the local CR spectrum is famously exhibiting a power law of spectrum « E727 
up to an energy of ~3 PeV, where the slope becomes « E~°, this feature being usually referred 
to as the knee [4]. Let us here mention that some collaborations have also found that the knee 
of the light component (He+protons) might be located at lower energy around ~0.7 PeV [5,6], 
which could mean that the proton knee is below 1 PeV. 

The remarkable, and almost featureless power law (although to be precise several instru- 
ments have in fact revealed some deviations in the powerlaw [7]), up to the knee, is seen as 
an indication that the sources of Galactic CRs must be accelerating particle up to the knee, 
and therefore, be pevatrons. Hence, the wild hunt for Galactic pevatrons, with a sometimes 
oversimplified argument that “finding the pevatrons” would mean “unveiling the origin of 
CRs up to the knee”, or “unveileing the origin of Galactic CRs”. Let us here insist on the fact 
that for the source(s) of Galactic CRs, being (a) pevatron(s) is a necessary condition but is not 
sufficient. It is for instance clear now that the Crab nebula is an astrophysical site where the 
acceleration of PeV particle (at least electrons) is taking place, and is however, not currently 
accepted as a typical source of the bulk of Galactic CRs (protons). 


1.2. 100 TeV Gamma Rays 


The best way to probe particle acceleration up to the PeV range, is (so far) through 
gamma-ray observations, in the very-high-energy range (VHE, 100 GeV to 100 TeV), and 
ultrahigh-energy range (UHE, above 100 TeV). Indeed, radiations in the VHE and UHE range 
directly testify to the abundant presence of nonthermal particles. Three mechanisms can 
essentially produce gamma rays in this range. First, a hadronic mechanism: the creation 
and decay of neutral pions through the interactions of accelerated hadrons with nuclei of 
the ISM. Second, a leptonic mechanism: the inverse Compton scattering (ICS) of accelerated 
electrons on soft photons (CMB, optical, infrared). Third, another leptonic mechanism: the 
Bremsstrahlung emission of accelerated electrons that can also lead to the production of 
gamma rays in this range [8,9]. 

In the ~TeV domain, it is not easy to differentiate between these different potential 
origins, and to know if parent hadrons or leptons are responsible for the production of gamma 
rays. The interpretation of the origin of observed TeV gamma rays is often subject to fervent 
discussions, and in many cases, the gamma-ray emission can be explained by one mechanism 
or another. See e.g., the numerous discussions on many well-studied astrophysical objects, 
for which the GeV and TeV gamma-ray emission can often be accounted for by a leptonic 
ICS origin, a hadronic origin, or a mixture of both. The case of the superova remnant (SNR) 
RXJ1713-3946 is a stereotypical example of the copious discussions on the interpretation of the 
gamma-ray emission. [10-22] 

It has often been claimed that above 250 TeV, the interpretation of detected gamma rays 
becomes less ambiguous, since the inverse Compton scattering of nonthermal electrons be- 
comes dramatically inefficient in producing gamma rays, due to the Klein—Nishina suppresion: 
at these energies, gamma rays are thus expected to testify that the acceleration of hadrons is 
taking place. Moreover, the VHE gamma rays of hadronic origin mimic the parent proton 
distribution, with a typical energy of gamma rays scaling as E} ~ 1/10E,. Therefore, this has 
often been condensed into “the detection of ~100 TeV gamma rays is be a direct indication of 
acceleration of ~PeV CR protons”. 
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Still, this affirmation needs to be tempered. The distinction between hadronic and 
leptonic origin of gamma rays is far from being straightforward even in the ~100 TeV range. 
Although the Klein—Nishina suppression considerably reduces the production of gamma rays 
from ICS in this range, leading to an exponential suppression above 250 TeV, gamma rays can 
still be produced by electrons through ICS. A conspicuous amount of VHE electrons can thus 
produce a gamma-ray signal in the 100 TeV region. Such a gamma-ray signal can especially 
be important in the context of instruments operating with an ever increasing sensitivity 
in the 250 TeV range, hence capable of detecting gamma rays photon even from steep or 
exponentially suppressed spectra [23-25]. Moreover, the synchrotron emission of accelerated 
electrons can also efficiently produce gamma rays. 

A prime example of a gamma-ray emission interpreted as the result of accelerated elec- 
trons is the Crab Nebula, the brightest source of the TeV sky, and thus one of the best studied 
objects [26-28]. More generally, several theoretical works have illustrated that electrons accel- 
erated above 2100 TeV could lead, through inverse Compton scattering or bremsstrahlung, 
to gamma rays in the 2100 TeV range (see e.g., [9] and references therein), and can produce 
hard gamma-ray spectra extending up to the = 100 TeV, making it especially strenuous to dif- 
ferentiate between hadronic and leptonic acceleration mechanisms [29]. Let us additionnally 
mention that the detection of neutrinos would be an effective and elegant way to discriminate 
between the hadronic and leptonic mechanisms since neutrinos can only be produced through 
hadronic interactions [30]. 

The detection and study of electron accelerators is of course of great importance, but in 
the search for the origin of CRs, is not enough. First, because CRs are mostly protons, and 
we are thus looking for proton sources: an efficient electron accelerator is not necesseraliy an 
efficient proton accelerator. The case of SNRs is instructive, since in spite of extensive studies 
the content of electrons and protons accelerated at SNRs is still poorly understood. Second, 
because unlike protons, electrons suffer important losses while transported in the Galaxy. 
Therefore, the problem of the origin of CR electrons is more local than the one of CR protons 
and somewhat disconnected [31-34]. 

So far, on the observational side, all major observatories relying on Imaging Atmospheric 
Cherenkov Telescopes (IACTs), such as VERITAS [35], MAGIC [36], or H.E.S.S. [37], and 
predecessors, such as HEGRA [38], have been probing VHE domain up to a few tens of TeV. 
Therefore, the search and identification of pevatrons has been a remarkably challenging task, 
but nonetheless successful since H.E.S.S. was capable of claiming the detection of the first 
pevatron in the Galactic center region [39]. Other observatories, using different techniques, 
relying on shower front detectors (or “air shower detectors”), such as HAWC [40], Tibet 
ASy [41] and LHAASO [28], operating up to and above the ~100 TeV gamma-ray domain 
also reported on the detection of several Galactic pevatrons. All these observations indicated 
that most of the pevatron candidates seem to not be associated to SNRs. 

In this short review, we intend to especially discuss the case of SNRs, since they have 
been, for long, at the center of the debate on the origin of Galactic CRs, and that the question of 
whether they can accelerate PeV particles is still open. We then briefly discuss other pevatron 
candidates and important open questions in the coming years. For pedagogical introduction 
to the general physics of CRs, we refer the reader to reference writings of the field [1,42,43]. 
For updated reviews discussing recent results and advances on the problem of the origin of 
CRs, we refer the reader to Blasi [44], Gabici et al. [45]. 


2. Supernova Remnants 


In the search for the origin of Galactic CRs, supernova remnants (SNRs), the spherical 
shock waves expanding in the ISM after the explosion of massive stars, have become over 
the years the most famous candidates. Several reviews detail the reasons of the success of the 
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SNR paradigm [46-48]. In short, the SNR hypothesis is supported by at least the following 
strong arguments: 


1. 


Energetic considerations: a simple back-of-the-envelope calculation illustrates that con- 
verting a fraction (~1%-10%) of the total explosion energy of SNe into CRs can account 
the measured energy density of CRs at the Earth. 

Spectrum of accelerated particles: at SNRs, which are strong shock waves (i.e., of 
compression factor 24), a first order Fermi mechanism, referred to as Diffusive shock 
acceleration (DSA) [49-52] can explain the efficient acceleration of particles with a power 
law in momentum «œ p~“, with a ~ 4 for a strong shock (which corresponds, at high 
energy to E~“**), somewhat compatible with the Local CR spectrum measured at the 
Earth. 

Magnetic field amplification: several mechanisms have been shown to be able to amplify 
the magnetic field at SNRs shock fronts, therefore helping to the reach the PeV range. 
Indeed, the Hillas criterion, where the Larmor radius of particles rz is equated to the 
typical size of the accelerator (the SNR radius r,,) gives that the maximum energy of 
particles is typically at most: 


(Tsh Ush B 
RR (=) (m Ez) (Z) (1) 


with usn the SNR shock velocity, and B the magnetic field. Therefore, in order to attain 
the PeV range, for a typical SNRs of a few pc, expanding at a few 1000 km/s, values 
of at least > 10? uG are needed. More precisely, there are now established theoretical 
results on the growth of instabilities at colisionless shocks, coupled to observational X-ray 
measurements at SNR shocks, that have reported values of magnetic fields of several 
hundreds of uG, i.e., several orders of magnitude above the values of the ISM [53]. 
Gamma rays: the detection of GeV to multi-TeV gamma rays from SNRs, and especially 
the capacity of current IACTs to resolve the shells corresponding to the spherical expand- 
ing shock waves, is an undeniable indication that efficient particle acceleration is taking 
place. 


In spite of these compelling arguments, the SNR hypothesis is however facing several 


difficulties, discussed in detail in several reviews [45,54]. Amongst them: 


1. 


Gamma-ray spectra: in the TeV range, the spectral index of all SNRs appears to be >2. 
This is problematic since the gamma-ray spectrum is expected to mimic the spectrum of 
particles accelerated at the SNR shock. Hence, this seems to indicate that the spectrum of 
particles accelerated at the shock is steeper than what is expected in the test particle limit 
œ E~?. Although the authors of various works have studied how efficient CRs at shocks 
can modified the produced spectrum, though nonlinear effects [55-60], the question of 
the slope of particles accelerated at SNR shocks is still open. 

The spectrum of particles injected in the ISM: from B/C measurements, the effects of 
the transport in the Galaxy have been estimated, finding that it should affect the CR 
spectrum by « E~°3--0-6 (see e.g., [61]). In order to produce the spectrum of the local CR 
spectrum below the knee « E~”, this means that the CR sources must inject in the ISM 
a spectrum « E721..24 ie. steeper than what is predicted in the usual test particle limit 
(x E~?). The problem of the spectrum injected by SNRs in the ISM is still open [62,63]. 
The 27Ne/?°Ne abundance measured in CRs is a factor ~5 larger than the one in the solar 
wind [64,65]. This overabundance is not easily accounted for under the SNR hypothesis 
alone [66]. This has motivated the study of alternative scenarios, in which particles 
are at least in part accelerated in an environment enriched by stellar outflows (e.g., 
massive stars, stellar clusters, OB associations, superbubbles) and might explain the 
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overabundance of 7*Ne/?9Ne [67]. Other issues in the chemical composition of CRs have 
also been reported: for example, the enhanced *°Fe/*°Fe ratio, or the linear dependency 
of the B and Be abundance with metallicity in metal poor halo stars, that are not easily 
explained under the SNR hypothesis [54]. 

4. No pevatron: As mentioned in the introduction, the sources of CRs are expected to be 
pevatrons. This means that SNRs are expected, for at least a period of their evolution, 
to be pevatrons. Since the gamma rays from accelerated protons scale as E} « 1/10Ep, 
the detection of gamma rays in the ~100 TeV range probes the acceleration of PeV CR 
protons. So far, the observation of all SNR shells in the VHE domain have revealed cut 
offs indicating that PeV CRs are not efficiently produced. 


This last point is crucial, since of all the aspects mentioned above, the identification 
of a SNR pevatron is often seen to be the conclusive proof that SNRs are the sources of 
Galactic CRs. Hence, this has stimulated the investigation of alternative scenarios for the 
origin of CRs. At least two comments are in order here. First, stricto census, current IACTs 
are not sensitive enough to detect gamma rays in ~100 TeV range, but they still have been 
able to measure the gamma-ray spectrum in the TeV and tens of TeV region, and revealing 
exponential suppression that indicate that the flux in the 100 TeV range should be at least 
drastically suppressed. Second, recent results of several instruments, HAWC, TibetASy and 
LHAASO [28,68,69] have reported on the detection of ~100 TeV gamma rays from a region in 
which the SNR G106.3+2.7 is located. By design, the improved sensitivity in the 100 TeV range 
of these observatories, compared to typical IACTs, comes witha poorer angular resolution than 
typical IACTs. It is therefore difficult for these instruments to spatially resolve the source(s) 
of the 100 TeV emission and specify which astrophysical object(s) is (are) responsible for the 
acceleration of particles up to the PeV range [68]. In addition, the associated SNR, G106.34+2.7 
is middle-aged (~10 kyr) which makes it a priori a relatively poor pevatron candidate [70]. 
The role of the different objects located in this complex region have thus yet to be clarified. 

Two possible explanations for the fact that all known TeV shells seem to not be pevatrons 
are: either that SNRs are only pevatrons for a short period of time, a priori in the early stages 
of their evolution and that all studied SNRs are already too old to accelerate PeV particles, 
or/and, that only a small fraction of SNRs are pevatrons [71,72]. In both cases, this coud 
explain the limited chances to identify an active SNR pevatron, even in the future Galactic 
surveys. 


2.1. Acceleration up to the PeV Range 


One of the reasons why all best studied SNRs are not pevatrons might be because they 
are too old to still be pevatrons. Indeed, as a SNR blast wave expands in the ISM, it slows 
down. The maximum energy of accelerated particles is thus expected to decrease. Naively, if 
we estimate that the maximum energy can be estimated by equating the diffusion length of 
particles lą to a fraction x of the shock radius rsp: 


D(E 
la = P m aralt) 2) 


assuming a Bohm-like diffusion coefficient, for relativistic CR particles, D(E) = 4r,(E)c, with 
Bt, Tsh X t*, and ush & t*—1, this leads to: 


Emax (t) « B(£)ren(t)Usy(t) x #2416 ü 


thus Emax decreases in time provided that B > 2a — 1. For a typical SNR expanding in a 
uniform ISM, « = 4/7 [73,74] thus Emax decreases if B > 1/7. For a SNR expanding in a wind 
ar +, a =7/8, B > 3/4. 
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The nature of the mechanism driving the amplification of the magnetic field, needed 
to reach the PeV range is currently still a matter of debate. It has been shown that several 
mechanisms can theoretically lead to the growth of instabilities. Let us for example mention see 
also the short review [75]: 


1. Acoustic instabilities, where a gradient of the CR pressure VP in a fluid of density 
p leads to the acceleration of the fluid x VP/p. In an inhomogeneous medium, the 
inhomogeneities of density will in turn lead to inhomogeneities in the acceleration, 
therefore amplifying the starting density fluctuations [76-79]. 

2. Density fluctuations in turbulent magnetized plasma, where instabilities are amplified 
by the dynamics of the shock (and where CRs do not play any role in the growth of 
instabilities): the production of PeV CRs by this mechanism is however expected to be 
rather limited [80,81]. 

3. Resonant streaming instabilities: induced by accelerated CRs, the background plasma 
reacts to moving CRs and creates a current compensating the positive charge excess 
created by CRs. For particle of Larmor radius r,, this leads to the growth of resonant 
waves of wave number k ~ 1/r,. The resonant nature of this process typically leads to a 
saturation when the amplified field 6B becomes of the order of the pre-existing ordered 
magnetic field Bo. [82,83] 

4. | Nonresonant streaming instabilities (often referred to as Bell instability): as in the res- 
onant streaming instability, the current of CRs jcr accelerated at the shock cause the 
background to react with an oppositely directed current. This results in a force jer x B 
responsible for the growth of perturbations in the magnetic field. The growth is the 
fastest on scales much smaller than the Larmor radius of particles, hence the nonresonant 
label [84-87]. 


This so-called nonresonant Bell mechanism is especially important in the case of SNRs, 
because it can lead to magnetic field amplification with values significantly larger than the 
pre-existing magnetic field Bg. Moreover, the growth rate of the instabilities is typically faster 
than in all previously mentioned mechanisms, and is therefore thought to be the mechanism 
governing magnetic field amplification at SNRs. A pedagogical discussion on the growth 
of resonant and nonresonant streaming instabilities can be found in [44]. The growth of 
the magnetic field from nonresonant streaming instability is exponential until saturation is 
reached. The exact details on when saturation is reached is still a matter of discussion, but the 
typical amplified magnetic field 5B can be estimated equating the Larmor radius of particles 
to the typical plasma spatial displacement induced by the growth of instabilities: 


5B(t) © antre (4) 


where ¢cr is the CR efficiency (i.e., the fraction of the shock ram pressure converted into CRs 
at the shock through DSA), p the density upstream the shock front, and A = In(pmax/mc), 
and the slope of accelerated particles is assumed to be « p~*. Equation (5) leads to typical 


values: oe i 
ayo f SCR Osh (t) n \1/2 
SEU 2( 01 ) (sai km/s G ze. Ha (5) 


Thus, for shock velocities of a few thousand km/s, and density sufficiently high, which 
is typical for instance of the shock launched by the explosion of core-collapse supernovae 
exploding in the dense wind of a late sequence massive star, values of few hundred uG can be 
reached. 

The growth of nonresonant instabilities drives the amplification of the magnetic field 
and thus allows to reach higher values for the maximum energy of accelerated particles. 
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Refs. [85,86] discussed that this maximum energy can be estimated by considering that the 
growth of instabilities saturates when the amplification reaches a number N of e-folding. 
For the maximum growth rate Ymax, corresponding to wave numbers k, this condition reads: 
Ymaxt © N. The value of N is not well constrained, but arguments in favor of 3 < NV < 10 
have been made [85]. The corresponding maximum energy reads: 


— 3 
Emax (t) x Tein a 6) 


Expliciting Equation (7) with numerical values: 


nak SCR val SP t n 1/2 
Emax(t) = 1 ( 0.1 ) (sat E) (w cat) Gr = EEN g 


where N = 5. This shows that in order to reach the PeV range at times of the order of 
the century, relatively high shock velocities and density are required. However, a careful 
calculation including the precise quantities associated to typical SNRs shows that in most 
of the evolution of typical SNRs - in most of the free expansion phase, and in the adiabatic 
Sedov-Taylor phase - the PeV range cannot be attained. Indeed, in SNRs from thermonuclear 
supernovae, high shock velocities 25000 km/s can be observed but usually associated to an 
ISM of the order of ~1 cm~°. On the order hand, for SNRs from core-collpase SNe (CCSNe), 
which correspond to the bulk of SNe, the higher densities provided by the dense wind of the 
late sequence progenitor star can be as high as ~103-104 cm~? (decreasing with « r~?) but 
thereby substantially decreasing the shock velocity below [500 km/s. It therefore seems that 
only SNRs from peculiar SNe, sufficiently energetic and/or exploding in dense winds are 
capable of efficiently producing PeV particles [87]. Moreover, the duration of the pevatron 
phase of these SNRs appears to be rather limited. By estimating the total content of protons 
produced by typical SNRs and peculiar pevatron SNe, it has been possible to use the local CR 
proton spectrum to constrain the rate of Galactic SNR pevatron to a fraction of the Galactic 
SN rate <3-4% vsyn % 0.1/century [72]. Moreover, since these peculiar SNRs are typically 
pevatrons on timescales of K one century, this corresponds to one SNR pevatron active for the 
order of ~1 century every few ~10,000 years, drastically limiting the chances of catching an 
active SNR pevatron. In other words, it could be that SNRs are indeed producing the bulk of 
CRs, but that the short duration of the pevatron phase prevent us from seeing one in activity. 


2.2. Supernovae 


Discussions on magnetic field amplification, e.g., through the growth of nonresonant 
streaming instabilities, lead us to the idea that the maximum energy of accelerated particles 
scales in time as « tv3 (t) 2(t). This corresponds, for CCSNe in which the velocity is 
typically in the early stage (free expansion phase) vg, « t~!/8 in a wind of profile n œ r~? to 
Emax « t714. In the case of a thermonuclear SNe expanding in a uniform ISM, Vsh « 473/7 
in the free expansion phase, thus Emax « t~?/’. This indicates that the maximum energy 
is reached at the earliest time and supports the idea that if some SNRs are pevatrons, they 
most likely must be pevatron when they are the youngest. Since the duration of the potential 
pevatron phase is virtually unknown, some authors have investigated the possibility that 
PeV acceleration would especially be efficient during a few days/month/years after the 
explosion of SNe [88-99]. Consequently, such efficient particle acceleration could lead to a 
considerable emission of gamma rays from the GeV to the multi-TeV range, a priori mostly 
through hadronic interactions because of the high density of the circumstellar environment in 
which CCSNe explode. Of course, in our Galaxy, the problem of low number of events remain, 
with a SN rate typically inferred ~3/century, but it has been argued that the detection in the 
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gamma-ray domain of close-by extragalactic SNe could help to study the acceleration up the 
PeV range in the first days after the explosion of the SN (in other words, extremely young 
SNRs). 

The possibility for instruments to detect such SNe is essentially limited by a physical 
process that can degrade the gamma-ray signal: the two-photon interactions in which a 
gamma-ray photon interacts with a low energy photon from the SN photosphere, producing an 
electron-positron pair, can potentially degrade the gamma-ray by several orders of magnitude. 
In the case of luminous supernovae SN1993J, it has been shown that a detection of gamma 
rays from extragalactic energetic SNe within a few Mpc is possible with next-generation 
instruments [100]. 


2.3. Archeology of Pevatrons 


Whether it is SNe, or SNRs, or even other astrophysical objects, since the duration of 
pevatron phase is expected to be limited (month, years or few centuries), and that together 
with the low rate of these objects, we might be in a situation where the number of active 
Galactic proton pevatrons is low. See e.g., discussion above, on the number of SNR pevatrons 
from energetic SNe expected to be <0.1/century typically active for a duration of the order of 
~1 century. This would make it possible for us to be in a situation where SNRs do produce 
the bulk of Galactic CRs, but where we currently do not have any active proton SNR pevatron 
in the Galaxy. In such a situation, we need alternative ways to probe the past activity of 
pevatrons, and therefore attempt to do archaeology of pevatrons '. The signature of the 
extinct pevatron could for instance be found through deeper studies of SNR shocks (e.g., 
extrapolating backwards the shock dynamics), by probing the surrounding molecular clouds 
(MC) of the pevatron candidates, for instance through gamma rays observations of MCs in the 
vicinty of pevatron candidates [101-103], or through the study of the radiation of secondary 
particles produced around pevatrons [104]. 


3. The Detection of Galactic Pevatrons with Gamma Rays 
3.1. The Crab Nebula 


The Crab Nebula is probably one of the best studied objects of the Universe: its emission 
has been monitored over several decades and accross the entire accessible electromagnetic 
spectrum [105]. Extensive observations of the Crab Nebula have revealed variability in the 
gamma-ray domain [26,106,107]. These flares—typically lasting for several hours to few 
days—have often been interpreted as the result of efficient acceleration of electrons up to 
the PeV range, emitting gamma rays through synchrotron mechanism. Discussions on the 
acceleration mechanisms usually invoke magnetic reconnection around the Crab pulsar and 
shock acceleration at the pulsar-wind termination shocks [108-115]. Recently, the LHAASO 
instrument reported on the detection of gamma rays above =1 PeV, additional confirmation 
that an electron pevatron is located within the Crab Nebulae [28,116]. This means that, 
unexpectedly, the first detection of a pevatron corresponded to an electron—pevatron, a 
reminder that the finding of a pevatron does not immediately solve all problems of CR 
physics. 

Currently, the mechanism involved in electron acceleration at the Crab Nebula are still 
not well understood, which makes the case of the Crab Nebula even more remarkable. Several 
mechanisms have been proposed, such as for instance: (a) DSA at the termination shock 
formed between the pulsar relativistic wind and the circum-pulsar medium; (b) driven mag- 
netic reconnection of the alternating field lines compressed at the shock; (c) resonant cyclotron 
absorption in ion-doped plasma, in which electrons and positrons absorbe energy emitted by 
ions in the plasma (see e.g., for a short and pedagogical discussion [117]). Theoretical works, 


Universe 2021, 1, 0 


9 of 20 


through Magnetohydrodynamics (MHD) and particle-in-cell (PIC) simulations are helping to 
advance our understanding on particle acceleration around pulsars. 


3.2. The First Galactic Pevatron Detected through VHE Gamma Rays 


While the multiple observations in VHE range with IACTs of SNR shells were revealing 
that all known SNRs seem to not be pevatrons, observations with H.E.S.S. were objectifying 
that gamma rays up to ~100 TeV seem to originate from the Galactic center region around 
Sagittarius A*, a source called HESS J1745-290, thereby indicating the presence of a pevatron 
in the Galactic center. No SNR seem to be associated to this emission, making the first catch of 
a pevatron even more startling than expecting. As previously mentioned, the 100 TeV range 
is not directly accessible to H.E.S.S., but the observations were capable of revealing a hard 
spectrum with little suppression up to 10-20 TeV, strongly suggesting that 100 TeV gamma rays 
are produced [39]. This result as been confirmed by other instruments such as VERITAS [118]. 
Various interpretations of the gamma-ray emission of this pevatron candidate have been 
proposed, including e.g., supernovae and clusters of massive stars [94,119], millisecond 
pulsars [120-122], or particle acceleration from the supermassive black hole [123,124]. 


3.3. A Population of Galactic Pevatrons 


In the recent years, refined observations with various techniques, such as IACTs (e.g., 
H.E.S.S. [125], VERITAS [35], MAGIC [126]) (see e.g., for reviews [127—129]), shower front 
detector instruments (e.g., HAWC [130], Tibet AS [41], LHAASO [23]), and systematic survey 
of the sky, have helped discover several Galactic pevatron candidates. The HAWC observatory 
has for instance reported on the detection of at least ten sources with emission above 250 
TeV of likely Galactic origin [40,68,68]. The case of J2227+610 is especially exciting as a new 
potential Galactic pevatron [68], since it has also been detected with Tibet ASy [69], and its 
association with the SNR G106+2.7. Arguments in favor of the acceleration of PeV particles 
from the pulsar wind nebula VER J2227+608 through gamma-ray observations with Fermi- 
LAT had already been made [131], and make this case of special interest, although the situation 
remains quite unclear due to the complexity of the region [132,133]. 

Finally, with the detection of gamma rays of energy 2100 TeV by LHAASO from a dozen 
of sources that are likely located in the Galaxy, another important milestone has been reached: 
the first detection of a population of pevatrons [28]. So far, the spectral shape of most of these 
sources is still not well established, which makes it difficult to clearly understand whether 
the origin is leptonic or hadronic. Moreover, the source association remains a difficult task 
given the typical angular resolution of the order of ~0.3° at 100 TeV, but these detections are 
still an important step in the hunt for Galactic pevatrons. In Table 1 we give a quick overview 
of the known Galactic pevatron population as of May 2021, although this list is likely to be 
lengthened soon due to ongoing surveys. 


Table 1. List of known Galactic pevatrons as of May 2021. This list is likely to be lengthened soon due to active ongoing detection 


campaigns. 

Source Possible Association Reference 
HESS J1745-290 Sagittarius A* /Galactic center [39] 
Crab /LHAASO J0534+2202 PSR J0534+2200 [26,28,41,107] 
LHAASO J1825-1326/2HWC J1825-134 PSR J1826-1334/PSR J1826-1256 [28,134] 
LHAASO J1839-0545/2HWC 1837-065 PSR J1837-0034/PSRJ1838-0537 [28,40] 
LHAASO J1843-0338/2HWC J1844-032 SNR G.28.6-0.1 [28,40] 
LHAASO J1849-0003 PSR J1849-0001 /W43 [28] 
LHAASO J1908+0621/MGRO 1908+06/ SNR G40.5-0.5/PSR 1907+0602/PSR 1907+0631 [28,40] 
2HWC 1908+063 

LHAASO J1929+1745 PSR J1928+1746 /PSR1930+1852/SNR G54.1+0.3 [28] 
LHAASO J1956+2845 PSR J1958+2846 /SNR G66.0-0.0 [28] 
LHAASO J2018+3651 PSR J2021+3651/Sh 2-104 (HII/ YMC) [28] 
HWC J2019+368 [40] 
LHAASO J2032+4102/2HWC J2031+415 Cygnus OB2/PSR 2032+4127/SNR G79.8+1.2 [28,135] 
LHAASO J2108+5157 [28] 
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4. All Kinds of Pevatrons: Pulsars, Massive Stars, Stellar Clusters & Superbubbles 
4.1. Pulsars and Surroundings 


The Crab Nebula, as mentioned above, is a clearly accepted example of pevatron, where 
the gamma-ray flares reveal the acceleration of ultrarelativistic electrons, injected from the 
pulsar and are energized to at least the PeV range. Altough it is possible that protons are also 
accelerated, the Crab is a clear example of an electron pevatron. It illustrates the possibility 
for pulsars, and systems hosting pulsars to be efficient particle accelerator. The detection of 
several halos of TeV gamma rays, often referred to as TeV halos is connected to this issue of 
particle acceleration around pulsars [116,138,139]. 


4,2. Massive Stars, OB Associations and Stellar Clusters 


Aside from the fact that most SNRs seem to not be pevatrons, other issues remain: for 
example, the spectra of accelerated particles at SNR shocks, the spectra of particles injected in 
the ISM, or the electron-to-proton ratio accelerated and injected in the ISM. These tensions 
have motivated the search of potential alternative accelerators: massive stars, and their dense 
winds, leading to the formation of a strong shock at the termination of the shocks, were for 
instance proposed as Galactic CR factories [140-142]. Moreover, as massive stars tend to 
form in clusters, it has been proposed that the shock created by the collective effects of the 
winds of individual stars could accelerate protons above the PeV range, and contribute to the 
production of Galactic CRs [143,144]. For a review on particle acceleration from star forming 
regions we refer the reader to Bykov et al. [145]. 

In addition to theoretical arguments, the detection of gamma rays from star clusters has 
been a direct confirmation that efficient particle acceleration is taking place: Westerlund 1 
and Westerlund 2 detected by H.E.S.S. [146,147]. Moreover, spatially resolved observations 
of stellar clusters have revealed that the distribution of particles from massive stars seem 
to scale as the inverse of the distance away from the sources, which is an indication of 
steady production and injection of particles in the ISM. This is thus supporting the idea that 
these objects could be prominent sources of CRs [148-150]. Futhermore, the detection of a 
population of at least 12 sources emitting gamma rays above 2100 TeV by the new instrument 
LHAASO, testifies of the detection of at least 12 pevatrons of Galactic origin, at reported in 
Table 1. Although it is yet difficult to draw any firm conclusion on the association to known 
objects, most of these sources seem to not be associated to SNR shocks. On the other hand, 
some sources, as for example, J2032+4102 could potentially be associated to massive stars. 
Precise measurements of the differential gamma-ray spectra of these sources are needed to 
understand if the emission is due to the acceleration of protons or electrons above the PeV 
range, and deeper studies of the spatial origin of this emission will help understand which 
astrophysical object are responsible for the production of these particles. 


4.3. Superbubbles 


In the line of star clusters, superbubbles, the giant cavities of hot and turbulent plasma 
inflated by repeated neighboring SNe explosion, have also been proposed as potential pe- 
vatrons. Superbubbles have especially been mentioned in the context of the search for the 
origin of CRs because they could potentially account for unexplained trends in the CR com- 
position. Indeed, measurement of the CR spectrum have indicated that the abundance of 
volatile elements (e.g., N, Ne, Ar) seem to correlate with the atomic mass; a correlation that 
is not clearly found for refractory elements (e.g., Mg, Si, Fe) (see [54] for a review on the 
topic) . It has been shown that a mixture of material with primordial solar composition and 
material enriched by stellar outflows or ejecta could explain the observed trends [151]: such 
a mixed ISM can typically be found in superbubbles. Moreover, superbubbles might help 
overcome another identified problem: the problem of the abundance of Be and B. In metal 
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poor halo stars, Be and B have been found to increase linearly with metallicity, rather than 
quadratically, which is expected in a standard scenario in which CRs are accelerated in a 
nonenriched ISM [152]. Superbubbles might help explain this linear dependency. This makes 
them interesting candidates in the CR origin problem, and as any good CR factory candidate, 
the question of particle acceleration up to the PeV range is central [153-157]. Detection of 
gamma rays from the Cygnus X superbubble [158] supports the need for further investigations 
with next-generation instruments optimized in the 100 TeV range. Finally, the detection of 
high-energy neutrinos from the Cygnus Cocoon region, in association with gamma rays, 
comes as direct indication that efficient acceleration of hadrons is taking place, and has been 
claimed to be a clear evidence that a proton pevatron is located in the Cygnus region [159]. 


4.4, Other Candidates 


In the list of potential proton pevatrons, and thus possible contributors to CR proton spec- 
trum, we must also mention other astrophysical objects, such as, for example, low-luminosity 
black holes binaries and their magnetically arrested disks [160], or direct acceleration of pro- 
tons in the magnetospheres of pulsar [161]. Although we only briefly mention these recently 
proposed scenarios here, they clearly deserve full attention and must be kept in mind in the 
hunt for pevatrons. 


4.5. Superpevatrons 


The improved performance of instruments have made possible the detection of gamma 
rays of energies even above 100 TeV. For instance, HAWC has reported on the detection of 
gamma rays around ~200 TeV from J1825-134 [68] and TibetASy on the detection of gamma 
rays 2100 TeV potentially associated to the region around SNR G106+2.7 or PSR J2229+6114, 
as previously discussed [69], and the LHAASO collaboration has reported on the detection of 
gamma rays of 21 PeV around the source J2032+4102, or from the Crab Nebula [116]. Should 
these gamma rays be due to hadronic interactions, this would indicate the acceleration of 210 
PeV CR hadrons, and thus point to Galactic superpevatrons. These superpevatrons might 
thus play a role in shaping the CR spectrum above the knee, and open many questions on the 
origin of CRs above the knee, and on the transition between CRs of Galactic and extragalactic 
origin. 


5. Gamma Rays: Limitations and Hopes for the Future 


Gamma-ray instruments have now clearly demonstrated their potential in the search 
for pevatrons. It remains that several or most of the aspects of particle acceleration up to 
the PeV range have yet failed to be understood through observations in the gamma-ray 
range: the astrophysical site(s), the physical mechanism(s), the duration of the pevatron 
phase(s). In order to move forward on these issues, it is essential to progress on the spatial, 
time and spectral resolution. This statement is quite general in Astronomy, since improving 
performance is always a strong motivation behind new generations of instruments at any 
wavelength: in the case of gamma-ray instruments and the search for pevatrons, it is now 
essential in VHE and UHE domain. It is worth mentioning typical orders of magnitude that 
are determining in the context of the hunt for pevatrons. First, the time evolution of the 
sensitivity of gamma-ray instruments makes it usual to target sources for at least several tens 
of minutes, and usual to need ~tens of hours, which naturally limits the amount of observing 
time available per source. In turns, this often prevents multiplying observations of a target of 
interest on monthly scales, and thus makes it challenging to probe the time dependence of 
the gamma-ray emission. Second, the angular resolution, even in the most favorable case of 
next generation IACTs (such as CTA) [24], is at best of the order of a few arcmin at ~100 TeV 
and ~10 arcmin at 100 GeV, which limits the possibility of probing extended emission even 
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around Galactic accelerators. This is problematic, for instance to address the question of the 
escape of particles from their accelerators, or to study complex regions around pulsars (e.g., 
TeV halos) [162]. 

Finally, the spectral perfomance. Again, even in the favorable case of [ACTs, obtaining 
a differential spectrum from a point source requires a significant amount of observing time 
(typically several hours), and the obtained spectrum is usually composed of ~10-20 bins per 
decade at ~1 TeV, with error bars increasing at high energy with the decreasing fluxes. A least 
two comments are in order. In the context of the search for pevatrons through gamma-ray 
observations, it is now clear that the shape of the spectrum above 50 TeV is crucial (and 
especially above 100 TeV), if only to discriminate between the potential leptonic and hadronic 
origin. We have seen that with the remarkable performance of air shower instruments (such 
as LHAASO), integrated spectra can reveal the presence of pevatrons, but the characterization 
of their spectra is also essential: spectral index, existence of an exponential suppression, 
gradual change of slope. These details imprinted in the gamma-ray spectrum can unveil 
precious information on the parent particle content. For instance, a sharp exponential cutoff 
in the 50 TeV range is expected in the case of parent electrons, due to the Klein—Nishina 
suppression, whereas a more gentle steepening could sign the presence of parent protons, 
or of synchrotron from electrons, as for example with the Crab gamma-ray spectrum and its 
gradual steepening over three energy decades [116]. Moreover, differences in the initial spectra 
of parent protons and electrons could also imprint the gamma-ray spectra. For example, at 
SNR shocks, it has been shown that the high-energy suppresion for protons is often well 
approximated by an exponential x exp(—p/ pmax) whereas for electrons, if the diffusion is 
Bohm-like and the maximum energy loss limited, the spectrum is expected to be suppressed 
as x exp(—(p/P%nax)”) [163,164]. 

Futhermore, it has been shown that the shape of the suppression at VHE energy can also 
vary for different accelerators, even considering only the acceleration of protons. This is for 
instance the case of stellar clusters. As detailled in Morlino et al. [144], which proposed a 
semianalytical description of the acceleration of particles at the wind termination shock of 
stellar clusters, the geometry, the diffusion regime both inside the stellar cluster wind and 
outside the termination shock can shape the suppression at the highest energies. In principle, 
the study of the gamma-ray spectral index at the highest energy could therefore offer precious 
indications on physical parameters at accelerators of interest, although, of course, such study 
will remain extremely challenging even for next generation IACTs. 

Eventually, let us observe that gamma-ray astronomers have taken the habit of trying 
to fit differential spectra with power laws, broken power laws, or power law exponentially 
suppressed, or other usual simple functions. This approach is physically reasonable, effective, 
and useful, but it can sometimes become decorrelated from any underlying assumption on 
the physical mechanism involved. Moreover, in many cases, the quality of the data after 
high-level analysis is still limited and allows for some freedom in the fitting. In some cases, 
one can then wonder if a given function providing a decent fit to a set of data, is really more 
motivated than another, and it has to be kept in mind that the fitting with usual function might 
be artificial and lead to loss of information. As mentioned above in the case of stellar cluster, 
at VHE, deviations from a clean exponential cutoff could sign different physical aspects, and 
it is especially important to keep this point in mind for future investigations of the gamma-ray 
spectra in the 250-100 TeV range. Indeed a wide variety of spectra could be revealed in this 
range, from steep to exponentially suppressed, or modified exponentials, including gradual 
changes of index or even bumps and features, holding valuable information. 

In addition to the gamma-ray domain, it is essential to mention observations at other 
wavelengths useful in the pevatron search, already down to the radio domain, where for 
synchrotron emission of accelerated electrons or of secondary particles can help probe the 
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processes at stake [71,104,165,166]. Multimessenger observations (especially neutrinos) are 
also encouraging [167-169] (see e.g. [30,170,171] for reviews). 


6. Open Questions and the Hunt in the Coming Years 


The recent detection of several Galactic pevatrons with various instruments operating 
the gamma-ray domain is gripping and might help us close in on the origin of Galactic CRs. 
Although it has often been claimed that the detection of 100 TeV gamma rays would be the 
unequivocal proof of the acceleration of PeV protons, and thus almost straightforwardly 
uncover the sources of Galactic CRs, it is now clear that things are a bit more abstruse. The 
detection of Galactic pevatrons has yet maybe brought more questions than answers. Among 
these questions: 


1. Is the gamma-ray emission from Galactic pevatrons due to the acceleration of electrons 
or protons? 

2. | What mechanism(s) drive(s) the acceleration of particle up to the PeV range? 

3. | What astrophysical object(s) is (are) responsible for the acceleration of PeV protons (and 
thus likely candidate for the origin of CRs)? 

4. How long do the pevatron phases last? 

5. When do PeV CRs escape their sources? What is the distribution of particles around the 
pevatrons? 

6. What is the population and distribution of Galactic proton and electron pevatrons? 

7. If the proton knee is, as suggested by some experimental results, below 1 PeV (~0.7 PeV), 
how does this affect the search for the sources of Galactic CRs? 

8. | Which objects are superpevatrons? How do they contribute to the CR spectrum? 


With the detection of several Galactic pevatrons, whose astrophysical nature and spectral 
properties have yet to be understood, it is now clear that finding 2100 TeV gamma rays, and 
catching a pevatron, is not enough to solve the problem of the origin of Galactic CRs. Sources 
of Galactic CRs must necessary be pevatrons, but this condition is not sufficient. In order to 
win the title of “main source(s) of Galactic CRs”, an astrophysical site will not only have to be 
a pevatron, but also to answer the different issues currently faced by SNRs such as: spectrum 
of particle injected in the ISM, chemical composition of CRs (77Ne/?2Ne ratio), fondamental 
physical process driving the acceleration of particles and magnetic field amplification, etc. 
Moreover, the role of the different detected pevatrons in the production of CRs and their 
contribution to the LIS will have to be understood. Fundamentally, one also has to prepared 
to the idea that not only one class (or “type”) of source(s) produce the bulk of Galactic CRs, 
but that the local CR spectrum might in fact the result of a complex intertwining of many 
sources, e.g., acceleration and injection at very high energy from some sources for a limited 
amount of time, reacceleration at other sources, etc., that would finally produce the measured 
local CR spectrum. The study of pevatrons will undoubtedly help to better understand the 
fundamental physical processes at stake (growth of instabilities, escape from sources, content 
injected in the ISM from sources), but will likely not be sufficient to answer all questions 
opened by CR physics. Complementary studies of the propagation of CRs in the Galaxy, or 
CR composition are vital to provide a complete understanding of the local CR spectrum. 

In the coming years, multiwavelength and multimessengers observations will also be 
essential to in order to progress on the question of origin of CRs. The gamma-ray range itself 
is expected to play a central role to probe the content of accelerated particles, and study the 
astrophysical objects involved. The instruments already operating in the VHE range, such 
as H.E.S.S., Veritas, MAGIC, LHAASO, HAWC, and the next-generation instruments, such 
as CTA [24], SWGO [25], with constantly improving spectral energy resolution and angular 
resolution, are needed to better understand the gamma-ray spectra in the ~100 TeV range, the 
emitting gamma-ray regions, and the distribution of particles around Galactic pevatrons. 
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